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Abstract 

■ We determine the probability that a random n by n symmetric matrix over {1,2,..., m} 
\ has determinant divisible by m. 

"S : 

1 Introduction 

S: 

'— — Let m be an integer. The m-rank of an integer matrix A is the greatest integer k such that ^4 
I \ has a by k submatrix (not necessarily contiguous) whose determinant is nonzero mod m, or 

£> ' if there is no such matrix. If m is a prime, the m-rank is equivalent to the usual rank over 

the field GF(m). In this paper we assume that the elements ay of A are chosen at random, 
^-j- | independently and uniformly, from Z m = {1,2,..., m}, subject to the condition that ay = a^j, 

l/l ' i.e., that ^4 is symmetric. For corresponding results without the symmetry constraint, see [1]. 

■^j- ■ Let P(n, m) denote the probability that a random n by n symmetric matrix A over Z m 

has m-rank n, and define Q(n,m) = 1 — P(n,m). Thus, Q(n,m) is the probability that 
| det(.A) = (mod m). As in Lemma 1.1 of [1], we have 

.^J ' Lemma 1.1 Suppose m = p^p^, 2 • • ~Pk k > where p±,P2, ■ ■ ■ ,Pk are distinctive primes. Then 

Q(n,m) = n Qtn.pf). 

8=1 

In view of Lemma 1.1, we restrict our attention to the case that m is a prime power, say 
m = p n '. It is useful to define q = 1/p. 

As in [1], our principal tool is Gaussian elimination, but in this case we have to use forms 
of Gaussian elimination which preserve symmetry. This is discussed in Section 2. Then, in 
Section 3, we use symmetric Gaussian elimination to show that P{n,p^) satisfies a five-term 
recurrence relation. In Section 4 we show that the five-term recurrence can be reduced to a 
three-term recurrence. Finally, in Section 5 we show that the three-term recurrence can be 
solved explicitly. 

The solution depends on the parity of n and fj,, and is well-known for fj, = 1, but appears to 
be new for fx > 1. To conclude Section 5, we deduce some inequalities from the explicit solution. 

An interesting problem is to determine the probability that a random n by n symmetric 
matrix A over Z m has given m-rank r, where r < n. The case fi — 1 has been solved by Carlitz 
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[2] (provided j) / 2), and the unsymmetric case has been considered in [1], but the general 
symmetric case remains open. 

Another open problem is to determine the probability that det(A) mod m takes a given 
(nonzero) value d. Small examples show that this probability depends on d. For example, if 
fi = 1, k = \n/2] , (d\p) is the Legendre symbol, 



s = 



0, if p = 2 or n is odd, 

(d|p)(-l) fc(p_1)/2 , otherwise, 



and IT is defined as in Section 5, then the probability is 



n 2fc (g) 



(l + sq k ). 



1-qJ U k (q 2 

2 Symmetric Gaussian Elimination 

Suppose that A,n,m = j/ 4 and P(n,m) are as in Section 1. Lemmas 2.1 and 2.2 describe 
symmetric versions of Gaussian elimination. 

Lemma 2.1 Suppose that n > 1 and an ^ (mod p). Then there is a matrix U such that 



( 



UAU 1 



an 



V o 



••• \ 



A' 



(mod p M ) 



(2.1) 



/ 



where A' is a random n — 1 by n — 1 symmetric matrix. 
Proof. Define 



U 



( 1 

A 2 1 

A 3 







1 J 



where 



o-nXj = —a\j (mod m) for j = 2, 3, . . . , n. 



(2.2) 



Observe that A 2 , . . . ,X n exist (since an ^ (mod p)) and (2.1) clearly holds. Also, A 1 depends 
linearly on the random symmetric matrix 



/ a>22 
V a n2 



a2n \ 



and thus is random over Z r 



□ 
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Lemma 2.2 Suppose that n > 2, an = (mod p), and a±2 ^ (mod p). Then there is a matrix 
V such that 



VAV T = 



an aw 





0-21 0,22 





A" 



(mod p M ) 



(2.3) 



where A" is a random n — 2 by n — 2 symmetric matrix. 
Proof. Define 

/ 1 

1 

A3 ^3 1 



V 



\ \ n \l n I J 



where 



an oi 2 \ ( Xj 

021 O22) 



\a 2j 



(mod m) 

for j = 3, ... , n. Observe that /Uj and Aj exist (j = 3, . . . , n) since 



det ( ° n ° 12 ) 7^ (mod p). 
\021 022/ 

It is easy to see that (2.3) holds. Also, A" depends linearly on the random symmetric matrix 

/ 033 • • • a 3n \ 



\a n 3 



Onn / 



and thus is random over Z r 



□ 



3 A five-term recurrence for P(n,p^) 

It is convenient to define 

P(0,p M ) = 1 for n > 
aild P(n,pV)=0 for^<0. 

Theorem 3.1 If n > 0,/x > 0, and boundary conditions are given by (3.1), then 

P(ra,p") = (1 - q)P{n - l,p") + q{\ - q^P^n - 2,p") 
+ q n {\ - q)P{n - l,^" 1 ) + q n+1 P(n,p^' 2 ). 



(3.1) 



(3.2) 
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Proof. Let A be a random symmetric n x n matrix over Z m , m = p^. The four terms on the 
right side of (3.2) arise from four mutually exclusive cases: 

1. an 7^ (mod p). 

2. an = (mod p) and some aij 7^ (mod p). 

3. aij = (mod p) for j = 1, . . . , n and an 7^ (mod p 2 ). 

4. ay = (mod p) for j = 1, . . . , n and an = (mod p 2 ). 

In case 1, which occurs with probability 1—q, we apply Lemma 2.1; since det(A) = an det(^4') 
(mod p M ) we have det(A) / (mod p^) iff det(A') 7^ (mod p' 4 ), which occurs with conditional 
probability P(n — l,p^). 

In case 2, which occurs with probability q(l — q 11 " 1 ), we can assume that ai2 7^ (mod p) by 
making a suitable permutation of rows and columns, if necessary. We then apply Lemma 2.2, 
obtaining det(A) = (a u a 2 2 - a\ 2 )det{A") (mod p»), so det(A) / (mod p^) iff det(A") / 
(mod p 11 ), which occurs with conditional probability P(n — 2,p M ). 

In case 3, which occurs with probability q n (l — q), we can reduce A to the form (2.1) 
because (2.2) is solvable. Thus det(A) = a u det(A') (mod p») and det(A) / (mod p^) iff 
det{A') / (mod p^ 1 ), which occurs with conditional probability P{n — l,^ -1 ). 

Finally, in case 4, which occurs with probability q n+1 , we can divide the first row and column 
of A by p, add random multiples of p^ 1 to elements in the first row (and to corresponding 
elements in the first column), add a random multiple of p^~ 2 to the (1,1) element, and obtain a 
new random symmetric matrix A such that 

det(A) = p 2 det(A) (mod p»), 

so det(^4) 7^ (mod p M ) with conditional probability P(n,p^~ 2 ). □ 



4 A three-term recurrence 

The five-term recurrence (3.2) with boundary conditions (3.1) can be used to calculate P(n,p^) 
in 0(n/j,) arithmetic operations. However, to obtain inequalities and asymptotic results it is 
useful to have an explicit solution. To obtain such a solution, we first reduce (3.2) to a three- 
term recurrence for P(n,p^) (n odd). 

Theorem 4.1 If p, > and boundary conditions are given by (3.1), then for odd n > 3, 

P(n,p") = (1 - q n )P(n - 2,j^) + q n P(n,p»- 2 ) (4.1) 

and for odd n > 1, 

P( «_ 1/) = «^!:yO. (42) 

Remarks. Equation (4.1) is a three-term recurrence from which P{n,p^) can be calculated for 
odd n. Equation (4.2) then determines P(n,p^) for even n. Equations (4.1) and (4.2) do not 
hold for all n > 3; for example, (4.2) fails if n is even and /j, is odd. In the unsymmetric case [1], 
(4.2) holds for both even and odd n. 
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Proof (of Theorem 4.1) Let 

oo 

V n = V n (x) = J2P(n,Pn^ (4-3) 

be a generating function for P(n,p M ). From the boundary conditions (3.1) we have 

Vo = x/(l-x). (4.4) 

Theorem 3.1 gives 

(1 _ q «+ia?)P n = (1 - q)(l + q n x)V n - 1 + q(l - q n - l )V n - 2 (4.5) 
for n > 1. Thus, for odd n = 2k + 1 > 1 we have 

(1 - q 2k+2 x 2 )V 2k+l = (1 - q)(l + tf tk+1 x)P 2 k + g(l - g 2fc )P 2fc -i (4.6) 
and for even n = 2k + 2 > 2 we have 

(1 - g 2fe+ V)P 2fe+2 = (1 - ?)(1 + q 2k+2 x)V 2k+1 + g(l - q 2k+1 )V 2k . (4.7) 

Both (4.6) and (4.7) hold for k > 0, and with the boundary condition (4.4) they define P n for 
all n > 0. 

Assume for the moment that (4.1) and (4.2) are correct. From (4.1) with n = 2k + 3 we have 

^ fe+ 3 = ^_ g2fc+ 3 x2 j (4.8) 

and from (4.2) with n = 2k + 1 we have 

^ = i _ g2fc+1 j (4.9) 

Clearly (4.8) and (4.9) for k > and (4.4) define "P n for all n > 0. Thus, it is sufficient to show 
that (4.8) and (4.9) for k > imply (4.6) and (4.7) for k > 0. 
From (4.8) and (4.9) we have 

which satisfy (4.6) and (4.7) with k = 0. Thus, we may assume A; > 1. Equation (4.8) with k 
replaced by k — 1 gives 

^-i = 1 _ q2k+1 J r 2k+1 . (4.i2) 

Substituting (4.9) and (4.12) in the right side of (4.8) and simplifying, we obtain (l-q 2k+2 x 2 )V 2k+ i, 
so (4.6) holds. Similarly, some algebra shows that (4.8) and (4.9) imply (4.7). 

Thus, the same generating function V n ,n > 0, is defined by (4.1) and (4.2) as by (4.5). □ 
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We can now give an explicit formula for the generating function V n . 
Theorem 4.2 If k > 0, then 

'l-q X 2 \ * / 



and 

(\- a 2k+1 x\ 

r 2k (x) = ( l _ q p 2fc+1 (x). (4.i4) 



, if fc >i. (5-1) 



Proof. Equation (4.13) follows by induction from (4.8), using (4.4) and (4.10). Equation 
(4.14) is just (4.9). □ 

5 An explicit solution and some bounds 

From Theorem 4.2 we can obtain an explicit solution for P(n,p fl ), and hence for Q(n,p^). Define 

n 

n„(g) = II(W) 
J'=l 

and 

f 1, if k = 0, 

Our explicit solution may be written in terms of T\ and T3: 
Theorem 5.1 // n > 1, // > l,k = [n/2\, and s = - 1)/2J i/ten 

P(n ' p,1) = (1 -"^(g 2 ) ( (1 " q2k+1)T ^ s) ~ 9/1(1 " 9 " )Tl(fc ' s) ) ' (5 ' 2) 

Proof. We may show by induction on k from (4.1) that 

P{2k + = (l^fiW) ( T s(*. s ) " s )) • ( 5 - 3 ) 

Considering the cases fi = 2s + 1 and /U = 2s + 2 separately, it follows from (4.2) that 

P(2k, P n = P(2k + l,p") + g 2 ^^^ Ti(k, s). (5.4) 

After some simplification we see that (5.2) holds both for n = 2k and n = 2k + 1. □ 

Equation (5.2) is inconvenient for numerical computation as P{n,p^) is close to 1 unless p^ 
is small. In order to deduce a convenient expression for Q(n,p^) = 1 — P(n,p^), we use the 
following identities for T\ and T3. 



6 



Lemma 5.1 If k > and s > then 



and 



Proof. (5.5) and (5.6) may be proved by induction on k. The proof is similar to the proof of 
Theorem 2.1 of [1], so details are omitted. □ 

If we substitute (5.5) into (5.2) the factor Tlk(q 2 ) /H 2 k(q) cancels. This gives a convenient 
explicit solution for Q(n,p^). 

Theorem 5.2 Ifn>l,n>l,k = \n/2\, and s = [(ji - 1)/2J then 

W) = -«">-« (5 . 7) 



where 



and 



R = g ^ £ (^(i - o - ^ +2 (i - ^)) ^ n^"^ ( 5 - 8 ) 



< P < g 3 ^ /2 . (5.9) 



Proof. Equations (5.7) and (5.8) follow from Theorem 5.1 and Lemma 5.1 after some simplifi- 
cation. Since jj, < 2s + 2 and n > 2k, we have q^(l — q n ) > q 2s+2 (l — q 2 i +1 ) for < j < k — 1, so 
P > (with equality only when k = 0). It is clear from (5.8) that R = 0(q fl+s+1 ) = 0(q 3 ^/ 2 ), 
and computation shows that R < q 3 ^/ 2 (the worst case is n large, fi = 2s+l large and q = 1/2). □ 

We now show that P{n,p^) has the expected monotonicity properties. 
Theorem 5.3 For n > and ji > 0, 

P(n + l,p M ) < P(n,p^) < P(n,^ +1 ). (5.10) 

Proof. The inequality P(n,p^) < P(n,p^ +l ) follows by induction onn + /i from the recurrence 
(3.2), since the coefficients on the right side of (3.2) are independent of fx. 

To prove P(n + l,p M ) < P(n,p M ) we consider several cases. If n is even the inequality follows 
from (5.4). If n = 2k + 1 is odd and fi = 2s + 2 is even then, from Theorem 5.1, 

P(2fc + l,p 2s+2 ) - P(2k + 2,p 2s+2 ) = g 2fc+2»+3 (1 + g) n 2fc+i(g) Tl(A . + l, s ) > 0. 

ii-kyq ) 

Finally, if n = 2k + 1 and /U = 2s + 1 are odd then, from Theorem 5.1, 
P{2k + l,p 2s+1 ) - P{2k + 2,p 2s+1 ) 

,.n fc+s (g 2 ' 
3 (<z 2 ) 

> o. 

□ 



^. + , + ^( ri(t + M) _ 9 ,5H 
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Corollary 5.1 lim n ^ 00 Q(n,p^) exists and lies in the interval q^/(l — q)\. Moreover, 

lim Pin.rr) = — — — 7T- > — ; ON — 

n^oo ^ > (l-g)noo( g 2 ) li^q 2 ) 

and 

where s = [(fj, — l)/2j . 

Proof. The limit exists by monotonicity in n, and the bounds follow from this monotonicity 
and Theorem 5.2. The explicit limits follow from Theorem 5.1 and Lemma 5.1 respectively. □ 
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